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Abstract. We report XMM-Newton observations of the eclipsing binary Algol A (B8V) and B (K2III). The 
XMM-Newton data cover the phase interval 0.35 - 0.58, i.e., specifically the time of optical secondary minimum, 
when the X-ray dark B-type star occults a major fraction of the X-ray bright K-type star. During the eclipse a 
flare was observed with complete light curve coverage. The decay part of the flare can be well described with an 
exponential decay law allowing a rectification of the light curve and a reconstruction of the flaring plasma region. 
The flare occurred near the limb of Algol B at a height of about 0.1 R* with plasma densities of a few times 
10^^ cm"'^ consistent with spectroscopic density estimates. No eclipse of the quiescent X-ray emission is observed 
leading us to the conclusion that the overall coronal filling factor of Algol B is small. 
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1. Introduction 

X-ray and radio observations have produced copious ev- 
idence for magnetic field related activity on esseritially 
all cool stars with outer convection zones fLinskv 1985: 
[Schmitt 1997; Schmitt & Licfkc 2002; Giidcl 2002). A 
comparison of solar activity indicators like X-ray lumi- 
nosity with those observed for other stars shows that the 
magnetic activity in the Sun is by comparison weak. Other 
stars can have X-ray luminosities exceeding that of the 
Sun by up to 4 orders of magnitude, yet this activity is 
usually interpreted by analogy with solar activity viewing 
the emission observed from stars as arising from scaled-up 
versions of the corresponding solar features. 

Because of the large distance of stars stellar coro- 
nae cannot be angularly resolved and therefore size in- 
formation is not available. The only exception comes from 
eclipsing binaries where size information can be obtained 
from light curve analysis. Particularly useful for coro- 
nal studies are those eclipsing binaries where one com- 
ponent can be considered essentially X-ray dark. Thus 
one has to look for eclipsing binaries whose primary com- 
ponent is of spectral type late B or early A. Only very 
few such systems suit able for X-ray analysis are known. 
In the system a CrB llSchmitt fc Kiirsteil l)l993l) ISchmitd 
l)l998|)lGiidel et all (|2003')) detected a total X-ray eclipse 
during secondary optical minimum, i.e., at the time when 
the X-ray dark A star is positioned wide in front of the X- 
ray bright B star, allowing them to reduce the size of the 
corona around a CrB B. Another such system is the pro- 



totypical eclipsing binary Algol, which consists of a B8V 
primary and a K2III secondary with an orbital period 
of 2.87 days. The relevant system parameters are listed 
in Tab. 1. Algol is among the apparently strongest coro- 
nal X-ray emitters and has been observed by all major 
X- ray satellites. Flares have bee n re ported, for example , 
bv lvan den Oord &: Me^ l)l989|) and lWhite et all l|l986>l 
who were the first to observe a complete (optical) sec- 
ondary minimum; from the absence of any eclipse at sec- 
ondary minimum, when the X-ray dark B8V star is in 
front of the K-type star, they argued for a corona with 
large scale heig ht around Algol B. On the other hand, 
ISchmitt fc Fav ata ( 1999) observed Algol for a whole bi- 
nary orbit using the Beppo-SAX satellite. They observed 
a giant X-ray flare with a decay time of approximately 
60 kiloseconds. At the time of optical secondary mini- 
mum a dip in the flare light curve was observed, which 
ISchmitt fc Favatal l)l999(l interpreted as an eclipse of the 
flaring plasma. From a detailed analysis of the Beppo- 
SAX X-ray light curve they constrained the size and lo- 
cation of the flaring region. Interestingly, the flare oc- 
curred above the south polar region of Algol B. Time 
resolved X -ray spectroscopy could also be obtained for 
this flare l)Fayata fc Schmitt! \l99^ . but from the low- 
resolution Beppo-SAX data no spectroscopic density di- 
agnostics was available. 

Here we report an observation of Algol with the EPIC 
camera and RGS spectrometer on board XMM-Newton. 
The observations lasted approximately 45 kiloseconds and 
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Table 1. System Parameters for Algol {— j3 Per) 



Parameter 


Primary 


Secondary 


System 


Mass (g) 


7.54 lO^"* 


1.64 10^^ 




Radius (cm) 


2.15 


2.29 10^' 




Spectral type 


B8V 


K2III 




Luminosity (erg/sec) 


5.95 10-''' 


2.44 10'''' 




Rotation period (days) 


2.8673 


2.8673 




Orbital period (days) 






2.8673 


Separation (cm) 






1.02 10''' 


Distance (pc) 






28 


Inclination (degrees) 






82.5 



10.0 



did cover the important phase interval between 0.35 and 
0.58, including optical secondary minimum. The XMM- 
Newton light curves have much better signal-to-noise than 
any previously reported X-ray light curve from Algol 
and have the specific advantage of a continuous cover- 
age, while, for example, the Beppo-SAX light curves were 
regularly interrupted by earth blocks. Further, the XMM- 
Newton observations were simultaneously accompanied by 
high-resolution spectroscopic observations with the XMM- 
Newton RGS, which in particular provide information on 
density-sensitive line ratios. 

2. Observations and Data Analysis 

The Algol system was observed by the XMM-Newton ob- 
servatory on February 12*'' 2002. The EPIC-PN detector 
operated with the thick filter inserted in the optical path 
in order to block optical light; data in full-frame window 
mode are available between 04:42:18 and 18:34:35 UTC. 
We extracted the EPIC PN lightcurve with standard tools 
provided by the Science Analysis System (SAS) software 
version 5.4. Although the source suffers from pile-up due 
to its substantial brightness (cf. Fig.^, we consider this 
pile-up effect as unimportant for the timing analysis pur- 
sued in this paper. In order to check for pile-up effects we 
compared the lightcurves applying different pile-up correc- 
tions and found no significant difference. We also checked 
the EPIC data for background proton flares, but none was 
found in the respective time interval. For our timing anal- 
ysis we extracted data from a circular region of 1 arcmin 
radius centered on the source while the background data 
was extracted from a source free region of the same size 
on the same CCD chip. A light curve was generated with 
the SAS tool evselect. Inspection of the background light 
curve showed it to be totally negligible compared to the 
observed Algol count rate; we therefore refrained from any 
further background subtraction. 

In Fig. E] we show the resulting XMM-Newton EPIC 
light curve vs. phase in four energy bands, i.e., below 
IkeV, in the band 1-2 keV, in the band 2-5 keV, and 
in the band between 5-10 keV. The phase refers to the 
phase of the primary minim um, which was c alculated 
from the ephemeras used by ISchmitt fc Favata y^22i)- 
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Fig. 1. EPIC PN light curves of Algol in the energy bands 
below 1 keV (squares), between 1-2 keV (dots), 2-5 keV 
(asteriks) and 5-10 keV (pluses). 

Interestingly, no change in the light curve is seen near 
the phase = 0.43, i.e., at the time of first contact, when 
the B type star starts occulting the K-type star. The light 
curves start to change at phase <j> — 0.48, when all light 
curves start increasing, but those at higher energy in a 
far more pronounced fashion. The light curves peak near 
the phase interval between 0.505 and 0.51, afterwards they 
start decreasing. At a phase (p — 0.5175 the decay of all 
light curves becomes very rapid, and a.t cj) = 0.525 the light 
curves in all four energy bands have reached pre-flare lev- 
els. Towards the end of the XMM-Newton observations at 
a phase of (f) — 0.555, the light curves in all energy bands 
start to increase again. 

An inspection of Fig.^shows that the light curves have 
the typical signatures of a stellar flare. There is a tendency 
for the higher energies to peak earlier in the flare also 
the flare spectrum signiflcantly hardens compared to the 
pre-fiare emission. For example, the flux between 5 and 
lOkeV increases almost by a factor of 10, while the flux 
below 1 keV increases only by 40% percent. However, the 
sharp decrease in the light curve between phases 0.5175 
and 0.525 appears to be very similar in all energy bands 
as well as the increase after phase 0.55. It is therefore 
extremely suggestive to interpret the dip in the flare light 
curve between phases 0.5175 and 0.55 as being caused by 
an eclipse of the flaring plasma by the X-ray dark B-type 
star, in complete analogy to the flare event reported by 
Schmitt and Favata (1999). However, in contrast to the 
Beppo-SAX observations of Schmitt and Favata (1999) 
the mid-eclipse was observed at phase = 0.50, while the 
mid-eclipse time of the XMM-Newton flare is at phase 
6 = 0.54. 




o 
01 

I 

c 

CL 

O 

Q_ 



0.35 0.40 0,45 0,50 0.55 

Phase 

Fig. 2. EPIC PN light curve of Algol in the total energy 
band; the solid line describes a simple exponential decay 
on top of a constant "background" emission; see text for 
details. 



In Fig. 12 we plot the light curve in the total energy 
band in a linear representation. The solid line shows an 
exponential decay on top of an assumed constant "back- 
ground" emission, i.e., a light curve with the functional 
form F(t) — Aq e~~, with Aq denoting the amplitude 
and r the decay time of the flare. It is apparent that the 
light curve after the flare peak at phase 0.51 can be well 
described by a such simple exponential decay. Exponential 
decays are very often observed for long duration solar and 
stellar flares, and seem to provide a convenient empirical 
form to describe flare decay light curves. For this flare 
on Algol we find fit parameters of Aq = 13.5cts/sec and 
T — 14.9 ksec. Note, that these values were determined "by 
eye" , rather than a formal fit procedure, and we refrain 
from giving a formal error analysis. We emphasize that the 
light curve data can be interpreted in such a fashion and 
that specifically the fiux level at the very end of the ob- 
servation is consistent with such an interpretation. Using 
these values we can rectify the flare decay light curve af- 
ter the flare peak by dividing the observed data points by 
the derived empirical flare decay law. The resulting light 
curve is shown in Fig. |3| By construction the light curve 
returns to the initial level after the end of the eclipse and 
the high symmetry of the flare eclipse light curve becomes 
apparent. 
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Fig. 3. Rectified EPIC PN fight curve of the flare dur- 
ing its exponential decay phase; solid line describes the 
modelled emission (see text for details). 



3. Reconstruction of the flaring region 

3.1. Methods 

The late start of the fiare eclipse indicates that the flare 
must have occurred near the limb of Algol B. From the 
system geometry it is also clear that the flare must have 
occurred on the "northern" hemisphere of Algol B, since 
otherwise the flare eclipse would have started prior to 
phase 0.5. In order to reconstruct the flaring volume we 
have to assume that the geometry of the flaring region did 
not substantially change during flare progress. Note, that 
the decrease in emission measure and hence flux has been 
taken into account by our rectification procedure, while 
changes in the geometry have not been (and cannot be) 
accounted for. Clearly, such an assumption is - strictly 
speaking - not valid. However, we point out that one may 
consider - at least approximately - the geometry of many 
solar flares as constant to first order. Without that as- 
sumption no further progress can be made in the stellar 
case. We next introduce a polar coordinate system with 
6 denoting the latitude measured from the north pole, 
0iong the longitudinal angle with the convention that at 
orbital phase (porh = 0.5 the longitudinal angle (/)iong — 
is oriented along the central meridian, and r denoting the 
height above the stellar surface. Obviously, for a limb flare 
there is a certain amount of degeneracy in 0iongi since 
one is looking almost tangentially through the corona. We 
therefore treat (/)iong as a fixed parameter which is adjusted 
to yield the most plausible reconstructed image. To span a 
volume we use a grid with iVi (=16) values in fj, — cos{9) 
and N2 (=13) values in , so that the whole grid has 
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■^tot = A^i X A^2 = 208 grid points denoted by some in- 
dex j = 1 . . . A'^tot- The light curve is given by Nph data 
points at phases 0^, denoted by an index fc = 1 . . . A^ph- 
With the system parameters as specified in Tab. ^ it is 
then possible to compute the visibility Vj^k of the volume 
element j at phase fc; Vj,k = 1, if the whole volume ele- 
ment is visible and contributes to the observed (rectified) 
count rate CRobs.fc at phase (j)k, and Vj,k — 0, if the whole 
volume element is invisible at phase (pk ■ Let Ij denote the 
intensity (measured in cts/sec), produced in the volume el- 
ement j = 1 . . . Ntot ; clearly, by construction Ij is constant 
and refers to the intensity at flare peak. The intensities Ij 
are related to the rectified observed count rates CRobs./c 
through the equations 



CRobs,fc = ^ ^j,klj, 



(1) 



and we are seeking that set oi Ij , j — 1 . . . iVtot , yielding 
a "best fit" to the observed rectified count rates CRobs,fc, 
fc = 1...7Vph, within their error ECRfc. Obviously, the 
solution of Eq.Q] represents a "standard" astronomical in- 
version problem. Specifically, Eq.^represents a discretized 
Fredholm integral equation, and such equations are well 
known to possess non-unique solutions. A good overview 
over inversion problems in astron omy and sol ution strate- 
gies for such probler ns is g i ven bvlLucvl (Il994') and the ref- 
erences given there. ILucvI l|l994l) specifically discusses the 
virtues of regularized and unregularized solution schemes 
and points out the importance of the constraint set by 
the positivity of the desired solutions. A rather trivial but 
nevertheless very important constraint is set by the con- 
dition 



/, > 0. 



(2) 



One inversion strategy involves minimizing the solution 
likelihood H defined as 



H 



E 

k=l 



CRcxp.fc ln(CRobs,fc): 



(3) 



where the quantities CRoxp,fcj k — 1 . . . iVph, are the "ex- 
pected" rectified count rates, which are calculated from a 
given current solution estimate Ij, j = 1 . . . Ntot through 
Eq.n This optimization problem can be iteratively solved 
by the algorithm ( Richardson.,,1972; ■Lucv..l994^) 



k=l CR„, 



CRo 



(4) 



starting with the initial estimate /j* — const. The problem 
with using the algorithm in Eq.^jis that there is no obvi- 
ous criterion when to stop the iteration. The solutio n like- 
lihood H continues to decrease on each iteration, but Lucvl 
l)l994j) shows inversions of problems with known solutions 
that the best fit reconstructed solutions tend to break up 
into delta functions, such that most of the "flux" ends up 
in a rather small number of bins. Thus the reconstruction 



of truly continuous d istribution functions becomes rather 
difficult. iLiicvl lll994l) shows that stopping the algorithm 
after a small number of iterations provides a good solution 
estimate and shows that a change in curvature in a HS 
plot does provide the derived stopping criterion. Here the 
entropy S is defined through 



Wtot 

i=i 



(5) 



ILucvI l)l994|) argues that such an unregularized reconstruc- 
tion is superior to regularized reconstructions because the 
choice of the value of the parameter "a" , which controls 
the strength of regularization, is arbitrary. We therefore 
decided to use Lucy's (1994) inversion scheme. 

3.2. Results 

Since the shapes of the light curves in different energy 
bands are very similar, we performed the reconstruction 
only for the total band. Next, in order to fulfill the positiv- 
ity constraint Eq. 13 we modelled the rectified lightcurve 
plus background, with an assumed constant background 
representing the emission from the rest of the star. The 
light curve used for reconstruction as well as the best fit 
model light curve resulting from an iteration of Eq. ^ is 
shown in Fig. 13 the reconstructed flare image is shown in 
a false-color representation in Fig. 0] In Fig. 0] the solid 
line shows the apparent limb of the K-type star, the dot- 
ted circles indicate heights in increments of 0.1 R*. The 
important feature of the reconstructed image shown in 
Fig. 21 and of all other reconstructed images considered 
by us is that the flaring region is rather well defined and 
confined. In particular, at the edge of the reconstruction 
volume, the reconstructed intensities become very small, 
indicating that the reconstructed size is not determined 
by the chosen reconstruction vo lume. 

As discussed extensively bv ILucvI l|l994j) . it is some- 
what difficult to decide when to st op the i terati on in Eq.^ 
We followed the recipe given by ILiicvI and com- 

puted for each iteration the likelihood H (cf. Eq. and 
solution entropy S (cf. Eq. In Fig. |Slwe plot S vs. H. 
There is no clearly identiflable point in the H-S diagram 
where the curvature changes; the solution chosen by us 
and shown in Fig. ^ is indicated by a large diamond. We 
also investigated the solution properties of further iter- 
ated solutions and found that those pixels with significant 
flux do not vary substantially in the course of further it- 
erations, while large flux changes are found only in pixels 
whose flux is so small that they do not substantially con- 
tribute to the overall emission. 



4. Discussion 

4.1. Flare energetics 

We first wish to consider the overall energetics of the 
flare. For this purpose we need to convert the observed 



J. H. M. M. Schmitt, J.-U. Ness, G. Franco; Algol limb flare 




20 



15 



10 




10 

Likelihood 



15 



20 



Fig. 4. Reconstructed spatial distribution of the flare re- 
gion on Algol B. All of the flaring plasma has been as- 
sumed to be located at a longitude of A = 70°. The circle 
represents the limb of the K-type star, units are in solar 
radius; dashed circles indicate heights in steps of 0.1 stel- 
lar radii. The system is shown at = 0.5556, when the 
longitude A = 70° is viewed at the limb. 

count rate into an energy flux. In order to convert from 
count rate to flux we considered various spectra of active 
stars (appropriately corrected for photon pile-up) and con- 
structed acceptable model fits to the EPIC data using a 
combination of individual temperature components. We 
found the energy conversion factor ECF = 5.4 10~^^erg 
cm^ ^ count ~^ to represent a reasonable choice of count 
rate to energy flux conversion (in the 0.2-10 keV en- 
ergy band) for active stars. With a distance of 28 pc to- 
wards Algol, we then compute a peak X-ray luminosity of 
6.7 X 10^° erg/sec for the flare. The total released X-ray 
energy Ex can then be computed from the observed de- 
cay time of 14.9 ksec as Etot ~ 3.6 x 10'^'^ erg; note, that 
this value only refers to the energy contained in the ex- 
ponential decay part of the light curve. Clearly, this flare 
is rather tiny with only a hundredth of the energy release 
observ ed in the giant flare observed bv lSchmitt fc Favatal 
l)l999(l . yet compared to solar flares, the flare is still about 
a thousand times more energetic than the strongest solar 
flares observed. As is clear from our reconstructed flare 
image, the actual size of the flaring region must have been 
quite small with a height not exceeding one tenth of a 
stellar radius. 

4.2. Size and emission measure of the flaring region 

The longitudinal size of our volume elements is - strictly 
speaking - undetermined. However, assuming that the ap- 



Fig. 5. Entropy (cf. Eq. EJ vs. hkelihood (cf. Eq. ISJ for 
the iterated solution. The chosen final solution is indicated 
by a large diamond. 

parent vertical size scale of the flare agrees with the longi- 
tudinal size scale, we deduce a longitudinal size of A 6°. 
Our light curve inversions do yield reasonable solutions 
also for other values than 0iong — 70°, however, the ac- 
tual value of 0iong cannot differ too much from that value. 
If 0iong was much larger, the K star would self-occult the 
flare, and also the rotation of the flare and the apparent 
motion of the B-type star would be in opposite directions. 
If '/'long was much smaller, the flare eclipse would start too 
early unless it is located far above the stellar surface. We 
therefore adopt - ad hoc - a longitudinal size of 10° for each 
of the volume elements, which then leads to a true size of 
Kicm = 1-0 X 10^°cm~'^ for all the elements. The recon- 
structed flare image (cf. Fig.^J) results in an effective count 
rate in each volume element. With Algol's distance of 
d=28 pc we can then compute the X-ray luminosity from 
each volume element through ix.oiom = CRcicmECF 47rd^. 
The optically thin cooling function for plasma tempera- 
tures at and above 10^ K is approximately given by P(T) 
~ ]^Q-24.75 rpo.25 gj.g(,j-|^3 gg(,-i^ j g^^ j-j- jg j.a^^;j^gj. tempera- 
ture insensitive. Adopting T 4.5 x 10^ K, which we de- 
termined fr om spec tral fits and also agrees with the values 
reported bv lYang~et a l. ( 2003|) . we can then compute the 
corresponding emission measure for each volume element 

through iJMelcm = Lxxlcm/P(T). 

4.3. Density of the flaring region 

Given Vdcm and EMdcm for each volume element, we can 
compute the plasma density ndom for each volume element 
using neiom = \/EMcicm/Vcicm- The resulting density fre- 
quency distribution in our volume elements is shown as a 
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Fig. 6. Cumulative density distribution histogram for re- 
constructed flare image shown in Fig. 01 From the recon- 
structed intensity in each volume element and its size, the 
emission measure and electron density can be calculated; 
shown is the density as a function of the integrated and 
normalized flare emission measure; see text for details. 

function of the total integrated and normalized emission 
measure in Fig. From Fig.|H|it obvious that the bulk of 
the flare emission comes from regions with high densities, 
more than three quarters of the emission measure resides 
at densities of at least 10^^ cm^^ and the peak densities 
reach approximately 2 x 10^^ cm^'^. With the simulta- 
neously taken XMM-Newton RGS data (with a slightly 
longer integration time of 52.66 ksec) we can check the 
densities derived from our light curve inversions. Doing 
this one must be aware that the bulk of the flaring plasma 
has a very high temperature (T w 45 MK), where essen- 
tially no Ovii emission is produced, and therefore Ovii 
emission actually provides no information on the main 
flare component. In Fig.[7|we plot the O vii Hc-like triplet 
obtained from the whole observation. The continuum is 
very high and the weak intercombination and forbidden 
lines of Ovii can hardly be detected above the contin- 
uum. We nevertheless extracted spectra during quiescent 
and flare phases, but we were unable to find significant 
differences in the f/i ratios. This indicates that either the 
sensitivity of our data is such that a density enhancement 
between flaring and quiescent emission is not detectable 
or that the f/i-ratio is actual ly dominated by th e radia- 
tion field of the B-type star ijNess et alJl2002alb|) . Given 
the geometrical configuration at secondary minimum, the 
observable corona on the K-type star is clearly fully im- 
mersed in the B star's radiation field. When accounting 
for the radiation field in the way described by iNess et all 
20n2lJ^ we estimate the density from our measured f/i 




21.5 



22.0 



22.5 



X (A) 



Fig. 7. O VII triplet obtained from the full XMM-Newton 
RGS data set; see text for details. 



ratio f/i = 0.86 ± 0.37 to log = 10.6 ± 1.42 (and 
otherwise log rig = H-l ± 0.23). We also note that the 
f/i-ratios derived for Algol from the Chandra MEG and 
LETGS and from the XMM-Newton RGS, i.e., from three 
independent data sets taken at three different times, all 
ag ree with each othe r (to within their respective errors; 
cf. iNess et al]l2003b|) . so that at present there is no evi- 
dence for phase and/or flare related changes in the Ovil 
f/i-ratio for Algol. Unfortunately, with the XMM-Newton 
RGS it is not possible to disentangle the contamination 
due to Fe lines in the Ne ix triplet (cf. iNess et a'l]l2003a(l 
and signal in the spectrum obtained from the second dis- 
persion order is too low to detect the Ne ix lines. Since the 
MEG and HEG f/i ratios derived from Chandra observa- 
tions at different times also do indi cate densities in excess 
of 10" cm-3 fcf. lNess et 'al]l2003bj) . we conclude that the 
densities derived from our light curve inversion are consis- 
tent with the available spectroscopic data on Algol. 



4.4. Cooling of the flare plasma 

An often encountered method to compute densities of a 
flaring plasma assumes radiative cooling to be the dom- 
inant cooling mechanism. Using the relation IPovle et alJ 
(1988) 



^e,rad 



3kT 

7p(t)' 



(6) 



where r denotes the observed decay time scale of the flare, 
we compute a density n^ rad ^ 9 x IQ^^ cm~^ leading to a 
volume of Vrad ^ 2 X 10'^^ cm'^, which is of a similar order 
as the volume derived in our reconstruction. We there- 
fore conclude that radiative cooling may in fact be the 
dominant cooling mechanism for this particular flare ob- 
served on Algol B , in co ntrast to the flare described by 
ISchmitt fe Favatal l|l999fl . which cannot have cooled by 
radiation alone. 
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4.5. Magnetic fields 

If we assume the post-flare material to be magnetically 
confined, we can compute the minimal field strength Bmin 
required for confinement from 

™" = 2n,kT. 



8tt 



(7) 



Using the values rie — 2 x 10^^ cm~'^ and T = 4.5 x 10 K 
as typical values for density and temperature, we find 
Bmin ~250 G. Assuming that the total released energy 
ultimately comes from reconnected magnetic energy, we 
can compute the typical reconnected, i.e., non-potential 
magnetic field strength from 



(8) 



d2 y 

reconnected 

Stt 

Using V — 5 X 10'^^ cm^, we find an annihilated field of 
about 430 G. Of course, since the observed released X-ray 
energy constitutes only some (small) fraction of the overall 
released energy, the latter and hence the total reconnected 
magnetic flux must be larger. If we assume - ad hoc - that 
only 10% of the released energy ends up as soft X-ray radi- 
ation, reconnected field strengths > 1000 G are required. 
We therefore conclude that coronal field strengths in ex- 
cess of 1000 G are present at least in the lower part of 
Algol B's corona. Thus a magnetic reconnection scenario 
for the Algol flare appears fully plausible. 



far more likely one. With the system parameters listed in 
Tab. m we compute that about 10% of the whole surface 
remain uneclipsed during the course of the XMM-Newton 
observations. Since the coronal densities of the quiescent 
emission is very likely (for a detailed discussion of dif- 
ferent density measuremen ts for A lgol and other 
coronal sources we refer to iNess et al...2003h) not 
too different from the density of the flaring region, which 
we know to be located at a height of ~ 0.1R+, we con- 
clude that the "quiescent" emission also comes from rather 
"compact" regions, which must, however, be located at 
latitudes 9 < 69° from the north pole. The "southern" 
eclipsed part of Algol B must have been devoid of any sub- 
stantial X-ray emission at the time of the XMM-Newton 
observations. Thus the filling factor of "quiescent" coronal 
X-ray emission of Algol B must be small, i.e, / < 0.1. The 
observed "quiescent" count rate of 25cts/sec corresponds 
to an emission measure of 1.3 x 10^^ cm~^, i.e., only about 
twice the emission measure of the flaring region. Assuming 
that the emission comes from the whole of the uneclipsed 
region, we compute for the product nlh = 2x 10 cm , 
where h denotes the average height of the corona. If that 
average height is the same as that of the flare, i.e, a tenth 
of a stellar radius, we compute a density of 3 x 10^° cm""^. 
Of course, the filling factor is likely to be smaller, and if 
the density is indeed 10^^ cm^"^, the overall filling factor 
on the K-type star is about one percent. 



4.6. Quiescent emission 

Finally, it is intersting to note, that at the time of first con- 
tact, which according to the parameters listed in Tab. ^ 
is expected to occur at phase (f> = 0.43 no change in the 
X-ray light curve is apparent. The light curve stays essen- 
tially constant until the fiare start a.t (p — 0.48. Eclipses of 
individual coronal background features might be hidden 
under the fiare-induced count rate increase, however, the 
smooth overall appearance of the fiare light curve provides 
little support for such an assumption. The mid-eclipse 
count rates in the energy bands 2-5 keV and 5-10 keV 
agree very well with those observed prior to flare and 
eclipse onset, so no eclipses need to be invoked. In the two 
lower energy bands the mid-eclipse count rates are reduced 
by about 10% compared to the pre-flare and pre-eclipse 
rates, on the other hand, the pre-flare and pre-eclipse rates 
are also somewhat modulated, so again the assumption 
of an eclipse does not appear to the urgently required. 
We therefore conclude that the "quiescent" emission from 
Algol B did not experience any substantial eclipse during 
our XMM observations and that the vast bulk of the "qui- 
escent" coronal emission remained unocculted during the 
eclipse, in contrast to the flare which was totally eclipsed. 

The absence of any clear eclipse of the "quiescent" 
emission can be explained cither by extended emission or 
by locating the emission within the unocculted north po- 
lar cap. As pointed out above, all of the available spectro- 
scopic density measurements of Algol indicate high den- 
sities, thus suggesting that the latter possibility is the 



5. Summary 

We have observed a flare on Algol B during secondary 
minimum with complete light curve coverage. The loca- 
tion of the flaring region will be well reconstructed with 
only some slight ambiguity as to the longitude of the flare 
occurrence. In particular, the height of the flare above the 
surface is about 0.1 and the flare did not occur in the 
polar regions. The plasma densities computed from the 
observed emission measure and the reconstructed volume 
agree with each other and are in the range of a few times 
10^^ cm^"^. The "quiescent" emission from Algol B was 
confined to a region at least 20° north of the equator, and 
only 10% of Algol B's surface were not occulted at any 
time during our XMM-Newton observations. If quiescent 
and flaring regions are located at the same heights, mini- 
mal densities of 3 x 10^*^ cm~'^ result; smaller filling factors 
result in correspondingly larger densities as suggested by 
spectroscopic measurements. Thus, the overall coronal fill- 
ing factor of Algol B appears to be quite small despite its 
enormous total X-ray luminosity. 
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